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We propose a highly efficient dual spin-valley filter in silicene, consisting of two distinct regions. In the first 
region, angular separation of the two valley spins in momentum-space is induced by a uniaxial strain, with 
further spin separation induced by an exchange field. The second region acts as an extractor of the 
requisite spin-valley current by means of localized fringe magnetic fields and gate modulation of the 
electrical potential. We demonstrated controllable and highly-efficient filtering (exceeding 90%) for all four 
spin-valley combinations based on realistic parameter values. We also discussed the feasibility of practical 
realization of the silicene-based spin-valley filter. 


The linear Dirac-like energy momentum dispersion of graphene has been instrumental for many 
spintronic and valleytronic applications. *“^^ The utility of the Dirac dispersion has motivated a simi¬ 
lar quest in silicene, which shares the same monolayer-honeycomb structure as graphene, but having 
heavier silicon atoms and consequently larger spin-orbit coupling.^^^^ 

In addition, silicene has the added attraction both theoretically^^ and experimentally^^ of real- and 
valley-spin dependence in its dispersion and transport properties, on account of its slightly buckled 
lattice. The inequivalent valleys K and K' at the corners of the reciprocal hexagonal lattice may be 
utilized in valleytronic applications. Recently, it was shown that application of strain which distorts 
the coupling strengths within the honeycomb lattice will also result in a valley-dependent gauge po¬ 
tential.^ **^ This valley-dependent effect of strain on the electrical properties of graphene has been 
investigated and proposed in several application such as the quantized valley hall effect, *^^ modula¬ 
tions of I-V characteristic of nanoribbons*^^ and valley filtering in graphene. *^^ In practice, the con¬ 
trollable strain can be generated by depositing onto stretchable substrates*^ *^^ and free suspension 
across trenches.**^ 

The dependence on both the valley-spin and real-spin degrees of freedom of electron transport in 
silicene suggests the possibility of inducing spin-valley polarized current in the material. Such spin 
and valley polarization of current has been achieved magneto-optically.*^^ An electrical method of 
inducing spin-valley polarization has been proposed by means of a Y-shaped spin-valley device.^**^ 
Similarly, a spin filter based on two dimensional U-shaped device^*^ and a three terminal Y-shaped 
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spin separator^^^ have been reported reeently. Another proposal involves the generation of spin-valley 
polarization by means of a slieene-based lateral resonant tunneling deviee,^^^ but this is effeetive only 
on normally ineident eleetrons. Finally, a teehnique^^^ was proposed to generate spin-valley eurrent in 
silieene by applieation of eleetrie field, whieh is however restrieted to eleetron energy eomparable to 
the relatively-weak spin orbit interaetion energy in silieene (of the order of a few meV). 

In this letter, we propose a double-barrier spin-valley filter based on strained silieene, whieh ean 
eleetrieally generate highly-effieient ( > 90%) spin-valley polarization of eurrent. The filter relies on 
material and operating parameters whieh are aeeessible experimentally. Unlike previous spin-valley 
polarization sehemes based on the spin-orbit eoupling effeet, our filter system is not restrieted to 
low eleetron energies of the order of the spin-orbit energy split Asoc (a few meV). The filtering 
prineiple is based on utilizing i) strain and exehange field to earve out distinet spin-valley transmission 
angular profiles in momentum-spaee, and ii) magneto-eleetrie fields in the seeond barrier to seleet the 
partieular transmission lobe eorresponding to the requisite spin-valley eombination. In this way, one 
aehieves a eontrollable means of generating any four of the spin-valley eurrent eombinations. 

We propose a two-barrier silieene filter shown in Fig. 1: Within the first barrier (0 < x < L), we 
apply a uniform uniaxial strain on the silieene lattiee, and an exehange field via adjaeent magnetie 
insulators on top and below the silieene layer, while in the seeond barrier {L + a < x < 2L + a), 
we apply a pair of d-funetion magnetie fields at the boundaries {x - L + a, and at x = 2L -i- a) via 
two ferromagnetie (FM) stripes, and modulate the eleetrieal potential via top and bottom eleetrostatie 
gates. In the model, these two regions take on different tasks: Basieally, the first region is utilized to 
ereate an angular separation (i.e., in 0-spaee, where (p is the angle of ineidenee) in the transmission 
of different spin and valley eurrents, while the seeond region is used as a barrier to sieve or filter the 
desired the eleetron eurrent of the desired eombination of spin and valley via the transverse Lorentz 
displaeement, and bloek all the others. For simplieity, these two regions are analyzed separately in the 
analytieal treatment. Subsequently, numerieal analysis are earried out on both regions in an integrated 
manner to verify the analytieal predietions. 

Under uniaxial strain in the x direetion and lattiee-dependent exehange field, the low-energy 
Hamiltonian ean be deseribed by 

Hrja-, = “ tlAskyjA + Wz^SOC + + l^cr- ( 1 ) 


In the above, rj = ±l denotes the valley index eorresponding to valleys K and K', respeetively, 
while (Tj, is the Pauli matrix eorresponding to real eleetron spin operator in the z-direetion. As shown 
in Eq. 1, the intrinsie spin-orbit eoupling energy term Asoc couples the real and valley spin degrees of 
freedom, assumes a relatively small value of 3.9 meV in silieene.^^^ The external spin-orbit eoupling 
(Rashba interaetion) is not taken into aeeount, as it is very small eomparing to the other effeets.^^^ 
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Fig. 1. Model of double-barrier, silicene based filter. Uniaxial strain is applied on the silicene region donated 
by yellow surface and hi,h 2 represent magnetic insulators (e.g., EuO) yielding to spin dependent energy gap 
in the region of (0 < ;<: < L). Second barrier consists of asymmetric FM stripes and electric potential barrier 
induced by top and bottom gates in the region of {L + a < x < 2L + a). Both barrier length L - 200 nm and the 
gap between the barriers a = 50 nm. 


The electrical potential in Eq. 1 is defined by /lo- = //q + where /iq is the shift in the potential 
due, e.g., to application of the gate voltage (in the first barrier /iq = 0 since there is no applied gate 
potential), while is the contribution due to the exchange field and is given by = 5 (^i + ^ 2 ). 
Here, /ti ,2 refers fo fhe exchange energy acting on fhe A and B sublaffices of silicene, which is induced 
by fhe fop and bottom magnetic insulator, respectively (Nofe fhaf due fo fhe buckling of fhe silicene 
layer, fhe A and B sublattice sifes have a vertical separation wifh respecf fo one anofher). Besides 
fhe spin orbif inferacfion energy, also fhere is an additional spin-dependenf effecl on fhe band gap of 
system due fo fhe exchange field which is given by A^ = \ {h\- / 12 ) cr^. However, fhis ferm is zero, 
since we have considered fhe case of parallel magnefizafion of bofh fhe top and bottom magnetic 
insulafors. Finally, application of sfrain in fhe Jc direction has opposite effecf on fhe local hopping 
energy, t ^ t + 6t {1 ^ 1.6 eV)^^^ for bonds in fhe fransverse y direction. This gives rise fo a sfrain 
gauge pofenfial of A 5 - 5t y which affecfs fhe ky componenf of fhe wave vecfor, fhus causing a 
fransverse deflecfion of fhe incidenf elecfrons. 


Lef us consider fhe expression for fhe election wave funclions in fhe fhree regions associafed wifh 
fhe firsf barrier: 


T/Cx) = e' 


ikxX 
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X < 0, 

0 < X < L, 

x> L. 


By subsfifufing fhe above ansafze info fhe Hamiltonian of Eq. 1 and evaluafing fhe corresponding 
eigenesfafes, one can evaluafe and Nj^a- in fhe above spinor expressions. These are explicifly 
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given by M„o- = ^+'?°'^soc ^ _ (£ /it^)+;?o-Asoc ^ p^j. simplicity cr - ±1 denote spin in the 

^ ^ ^ V(£-P.)2-('7^Asoc)^ ^ ^ ^ 

±z direetion. 

Next, we eonsider the wave veetors within and outside the first barrier: The x-eomponent of these 
are given by = ^|{E - - {qcrAsoc)^^^ and kjc = respeetively, where 

0(0) is the angle between the eleetron momentum and the x-axis within (outside) the barrier. By eon- 
sidering the eonservation of the momentum in the transverse y direetion in the unstrained (I and III) 
and strained regions (II), we obtain the relation: 


sin(0) Je^ - = sm{e)J{E -- {qcrAsocf + 


from whieh, one obtains the angle 6 = sin 


I I sin(0) sjE^-AgQi^+)]dt 


within the barrier. Finally, the valley 


\ Wt^tOM^rAsoc?/ 

and spin dependent transmission probability T" =\t^ ean be ealeulated by applying the wave funetion 
eontinuity relations at the boundaries. In Figs 2(a) and (b), it ean be elearly seen that the angular profile 
of fhe fransmission probabilify ehanges signifieanlly when fhe sfrain and spin-dependenf potential are 
varied. The dependenee of T on fhe sfrain field til ean be explained as follows: In fhe presenee of 
til, fhe eonservalion of fhe Iransverse momenlum al fhe inlerfaee belween fhe unslrained and slrained 
regions dielales lhaf E sin(0) = E sin(0) ± til, and Ihus, sin(0) = sin(0) + Therefore, fhe ehange in til 
direclly affeels fhe Iransmiffed angle of eleelrons. Besides, as slated earlier, Ihe momenlum qx is also a 
funetion of spin-dependenl potential Ha-. This resulls in furlher angular separation of Ihe fransmission 
probabilify T for dilferenl spin orienlalions. 

We found lhaf under Ihe optimal ehoiee of slrain and exehange slrenglh in Ihe firsl barrier, Ihe 
fransmission of parlieular spin-valley eombinafion ean be eompressed lo wilhin a lighf range of inei- 
denf angle 0. Moreover, fhe angles eorresponding lo perfeel fransmission (T" w 1) for all four differenl 
spin-valley eombinalions are well-separaled in fhe ineidenl angular spaee 0. This is a prerequisite for 
highly-effieienl filter operalion in fhe seeond barrier region. 

Lei us analyze fhe valley and spin separation wilhin fhe firsl barrier. Based on Eq. 2, and fhe fael 
lhaf perfeel fransmission in Klein lunneling oeeurs al normal ineidenee, i.e. sin(0) = 0, Ihe angles 
eorresponding lo perfeel fransmission for Ihe Iwo valleys are Ihen given by 0 = ? 7 sin“^ i i ■ 


Henee, one observes lhal Ihe separation of Ihe K and K' valley polarization in 0-spaee is indueed by 
Ihe slrain energy til. This simple analysis applies in Ihe limil of infinite barrier lenglh, i.e., L ^ oo. 
However, for a finite barrier Ihiekness, perfeel fransmission {T =1) oeeurs under Ihe resonanee 
eondilion: qxL = nn, where n = 0, ±1,.... For a sufiieienlly long deviee lenglh, exlra resonanee peaks 
oeeur eorresponding lo n > 2. In addition, Ihe dependenee of qx on Ihe real spin due lo Ihe spin- 
dependenl potential Her resulls in a furlher angular separation of Ihe perfeel fransmission direetion 
for opposite spins cr = ±1. These perfeel fransmission angles ean be derived analylieally for all four 
spin-valley eombinalions (cr - ± 1 , 7 / = ±l)by eonsidering bolh Ihe resonanee eondilion and Ihe 
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conservation of transverse momentum (Eq. 2). The analytical results agree with the location of the 
transmission peaks which are obtained numerically in Figs. 2(a) and (b). In the optimal case depicted 
in Fig. 2(a), the strain energy is chosen to be sufficiently large {6t - 43 meV) so that only the first 
resonant peak, i.e., qxL = ± 7 r occur, while the higher resonances correspond to angles greater than | 
and are thus totally reflected. Note that the strain energy cannot be excessively large, otherwise the 
transmission of even the first resonant peak will be suppressed. 

Having achieved a distinct angular separation of the transmission peaks for all four spin-valley 
polarizations by passage through the first (strained) barrier, we are in a position to effectively filter 
the desired spin-valley polarization by means of the second barrier. The spin and valley filtering are 
achieved within the second barrier by applying appropriate magnetic vector and electrical potentials. 
Delta d-function magnetic fringe fields B^{x) - BqIb{6{x) - 6(x- L)] are generated at the boundaries of 
the second barrier (here for simplicity we set x^O at the left boundary) by antisymmetric ferromagnetic 
stripes (see Fig. 1). These fields yield a square-hat vector (gauge) potential As = BolB[&{x)-Q{x-L)]y 
within the second barrier, where Ib = is the magnetic length. Essentially, this gauge potential 

induces a transverse deflection of the electrons, thus allowing a certain angular range in 0 -space to 
transmit through the second barrier. To filter out the requisite valley polarization, A b should match 
the magnitude of the strain gauge As - Physically, under this matching condition, the transverse de¬ 
flection from the strain A 5 and fringe field Ab cancels one another allowing electrons of the partic¬ 
ular valley index 77 to transmit through. From the matching condition, 77 ^ = sBoIb, one can obtain 
the required magnetic fringe field strength BqIb = ? 7 ^(^)^- In addition, to filter out the required 
real spin polarization, the electrical potentials of the two barrier regions should also be matched, i.e., 
A^o-.i = 7 Uo -,2 ^ crjUM,! = fJ- 0 , 2 , i-e-> the electrical potential 770 (due to the gate bias in the second barrier) 
should be equal to the (exchange induced) spin-dependent potential crfiM (in the first barrier). In this 
way, highly efficient filter operation can be achieved for the desired combination of spin (cr) and valley 
(? 7 )- 

So far, in our analytical treatment we have considered the two barrier regions separately. To verify 
our analytical prediction, we numerically calculate the overall transmission probability of whole 
system by considering the Dirac electron wave functions in all five regions and apply the proper 
matching conditions at the four boundaries. Subsequently, we calculate the conductance of spin-valley 
filter by integrating over the incident angles, i.e., 

Gfif = Go4cos(0)rf;fW, 


where Go = 


h^EpLy/hvF')' 


Finally, the spin-valley polarization (e.g., for valley K and spin t) is defined by 


p _ {GK1-GKl) + {GKt-GK'^)+{GKt-Gx'l) 

3{GKt+GKi+GK’t+GK'i) 
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PxT= 0-9412 Pk1= 0.9405 Pxt= 0.9412 PKi= 0.9405 


Fig. 2. The transmission probability of each spin and valley combination in terms of incident angle in the case 
of Ef - 25 meV, L - 200 nm, a - 50 nm, Asoc = 3.9 meV, a) polarization of different spin-valley in 0-space 
under the optimal configuration of strain and spin-dependent potential 5t - 43 meV, h\ - h 2 -2 meV as a result 
of the first barrier b) angular dependence of transmission probability under non-optimal configuration of strain 
and spin-dependent potential 6t - 16 meV, hi - h 2 - 2.5 meV. c), d), e) and f) show transmission probability 
of the system after passing through the second region. The magneto-electric potential in the second region is 
set at the optimal value so as to achieve an effective filter operation for all four spin-valley combinations: c) 
K-i, d) K^, e) Ki , and f) K.^. The magneto-electric configuration for the filter operations of the four spin-valley 
combinations are as follows: c) p = 2 meV and fitg = 43 meV, d) fi - -2 meV and 6b - -43 meV, e) yU = -2 
meV and Sts — 43 meV, f) — 2 meV and Sts — -43 meV. 
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Fig. 3. The analysis of the effect of varying strain field and spin-dependent potential due to exchange field 
/ii ,2 on the transmission Pk-\- In the following calculations, we assume the following parameter values: Ef - 25 
meV, L - 200 nm, a = 50 nm, Asoc = 3.9 meV. (a) The polarization Pk'i as a function of the strain field which 
is varied between 40 meV and 65 meV; (b) The dependence of the polarization Pk'\ on both the strain field and 
exchange field hx^', (c), (d), and (e) show the effect of changes in the strain field and the exchange field on the 
angular position of the peak transmission probability across the first region. The dashed black plot represents 
the angular position which matches the parameters (magnetic field and electrical potential) of the second barrier 
(filter) region. 


By assuming the optimal parameter values and matching conditions, we numerically obtain high 
valley-spin polarization > 0.94 for all four spin-valley combinations. Figs. 2(c) to (f) depicts the 
angular dependence of transmission probability of the whole system (transmission across 5 regions) 
under different matching configurations. In all four cases, we obtain transmission of only the desired 
spin-valley polarization while the other three polarizations are virtually filtered out. Additionally, one 
can control the transmitted spin-valley polarization merely by changing the sign of the magnetic and 
electrical potentials in the second barrier, i.e., by reversing gate-induced potential /iq ^ -po and the 
fringe magnetic field Our numerical calculafions assume a barrier widfh of L = 200 nm, and 

mafching magnefic field sfrengfh of Bq - 1 T and sfrain energy of df = 25 meV. These parameter values 
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are within reach experimentally, and thus the proposed spin-valley filter can be realized with current 
technology. By contrast, the spin-valley filtering solely by means of spin-orbit coupling effect^^^ yields 
a much greater overlap of the transmission profiles for fhe differenl valleys and spins, which franslafes 
fo a significanlly lower polarization values. Furfhermore, fhe valley-spin polarizafion exisfs only af 
low election energies, comparable fo fhe SOC energy of a few meV's. Numerically, if was found fhaf 
fhe SOC-induced valley and spin polarization affains a value of abouf 10% af election energy of 10 
meV, and vanishes when fhe energy is increased fo 20 meV. 

In Fig. 3(a), fhe valley-spin polarizafion Pk^ for fhe K valley and spin-up elections is ploffed over 
fhe range of sfrain field befween 40 meV and 65 meV. The plof shows a peak af 43 meV, fhe value 
which mafches fhe magnetic field sfrengfh in fhe second barrier (filler) region. For values of sfrain til 
slighlly above or below 43 meV, fhe polarization Pk^ decreases sharply due fo fhe mismalch of fhe 
vector potentials As and in fhe firsl and second regions. Inlereslingly, when bolh fhe sfrain til and 
fhe exchange field /ii 2 (and hence fhe spin-dependenl polenlial) are varied, we find lhal fhe reduction 
in fhe polarizafion Pk^ due fo increasing sfrain can be counler-acled by reducing fhe exchange field 
[see Fig. 3(b)]. This behavior can be understood by considering fhe effecl of Ihese Iwo parameters on 
fhe angular dependence of Iransmission probabilily. As shown in Fig. 3(c), when we sel til = 43 meV 
and h \2 - 2 meV, fhe peak of fhe Iransmission curve Tk^ (blue curve) exaclly coincides wilh fhe black 
dashed line which denoles fhe angular orienfalion which mafches fhe second filler region. However, a 
small change of til to 44 meV, resulfs in a significanl shifl in fhe Iransmission in 0-space [see Fig. 3(d)] 
away from fhe dashed line. This franslafes info a drop in fhe tiansmilled spin-valley of {K, t) across 
fhe filler region, and hence a lower Pk^. By modifying fhe exchange field to - 1 meV, fhe peak of 
fhe Tki curve can be made fo coincide wilh fhe dashed line [see Fig. 3(e)], and Ihus Pk-^ is reslored fo 
ifs peak value. 

We Iheorelically invesligale fhe real and valley spin polarized tiansporl in a double-barrier system 
based on stiained silicene, capable of achieving high real and valley spin polarizafion of currenf 
(exceeding 90%). The filtering principle is based on utilizing i) sfrain and exchange field in fhe 
firsl barrier to carve oul dislincl spin-valley Iransmission angular profiles in momenlum-space, and 
ii) magnelo-electiic fields in fhe second barrier to selecl fhe requisite spin-valley combination and 
filler fhe olhers. We demonstiale, bolh analytically and numerically, lhal almost pure spin-valley 
current can be realized based on geometric and material parameters that are within practical feasibility. 
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